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ABSTRACT: The Compact Neutron Spectrometer (CNS) is a Joint European Torus (JET) Enhance- 
ment Project, designed for fusion diagnostics in different plasma scenarios. The CNS is based 
on a liquid scintillator (BC501A) which allows good discrimination between neutron and gamma 
radiation. Neutron spectrometry with a BC501A spectrometer requires the use of a reliable, fully 
characterized detector. The determination of the response matrix was carried out at the Ion Accel- 
erator Facility (PIAF) of the Physikalisch-Technische Bundesanstalt (PTB). This facility provides 
several monoenergetic beams (2.5, 8, 10, 12 and 14 MeV) and a white field {E max ~ 17 MeV), which 
allows for a full characterization of the spectrometer in the region of interest (from ~ 1.5 MeV to 
~ 17 MeV). The energy of the incoming neutrons was determined by the time of flight method 
(TOF), with time resolution in the order of 1 ns. To check the response matrix, the measured 
pulse height spectra were unfolded with the code MAXED and the resulting energy distributions 
were compared with those obtained from TOF. The CNS project required modification of the PTB 
BC501A spectrometer design, to replace an analog data acquisition system (NIM modules) with 
a digital system developed by the Ente per le Nuove tecnologie, VEnergia e I'Ambiente (ENEA). 
Results for the new digital system were evaluated using new software developed specifically for 
this project. 



* Corresponding author. 
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1. Introduction 

Fusion facilities like the Joint European Torus (JET) require sophisticated diagnostic tools to effec- 
tively monitor the plasma inside the reactor torus. Numerous studies have been carried out using 
different types of apparatus, from simple optical devices to very sophisticated detector systems QU]. 
Still, to properly evaluate the reaction rate one needs to detect the products of the reaction; i.e., 
neutrons. Because neutrons interact little with matter, they are good candidates for properly char- 
acterizing the plasma. There are two neutron producing reactions which are of interest, the d+D 
and d+T reactions, which produce 2.45 and 14.1 MeV neutrons respectively. 

The Compact Neutron Spectrometer (CNS) project is a JET Enhancement Project for neutron 
fusion diagnostics in different plasma scenarios. It is based on a liquid scintillator (BC501 A) which 
allows for good discrimination of neutron (n) and gamma (7) radiation. The spectrometer is able 
to detect neutrons over an energy range that covers both of the neutron energies mentioned above. 
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To do neutron spectrometry with a BC501A spectrometer, it is necessary to have a reliable, fully 
characterized detector: it is crucial to have a good response matrix and, when operating under 
challenging circumstances such as the ones encountered at JET, the means to correct for the effect 
of any environmental conditions that cannot be effectively eliminated by shielding the detector. The 
Physikalisch-Technische Bundesanstalt (PTB) has many years of experience in the development, 
construction and use of liquid scintillation neutron detectors as neutron spectrometers The 
BC501A neutron spectrometer developed at the PTB has been shown to meet the requirements for 
fusion diagnostics [GL 01]. A system similar to the one described in this paper has been built by PTB 
and installed at ASDEX upgrade [Qj. The procedures for the determination of the neutron response 
matrix described here will be applied to the ASDEX upgrade neutron spectrometer as well. 

The experimental characterization of the CNS has been carried out at the PTB accelerator 
facility flfj. This experimental facility provides unique neutron radiation fields in the energy range 
from a few keV to 19 MeV. For the CNS project, a new method of data analysis had to be developed 
on account of the new digital acquisition system (DPSD) developed by the Ente per le Nuove 
tecnologie, VEnergia e VAmbiente (ENEA) which replaced the analog data acquisition system 
(NIM modules) that was previously in use Q70- The software necessary for the data analysis was 
developed at PTB. 

2. CNS Project 

The objectives of the CNS project can be summarized as follows: 

• Construction of the CNS at PTB and verification of the performance requirements. 

• Development of a computer code to carry out the analysis. The code was written in C++ and 
it makes use of the ROOT framework [ED- 

• Characterization of a BC501A fast-neutron detector with a digital acquisition system devel- 
oped by the ENEA group. 

• Determination of two response matrices (RM), an experimentally determined response ma- 
trix and a response matrix based on simulations adjusted to fit the data collected during the 
measurement campaign. 

• Evaluation of the resolution that is achievable after unfolding the measurements with the 
MAXED code, and comparison with Time-Of-Flight (TOF) measurements. 

Using unfolding techniques, the CNS should be able to achieve a resolution (AE/E) of 6% for the 
2.5 MeV peak and 4% for the 14 MeV neutron peak. It should also be able to handle the high count 
rates that are expected at JET [0], 

3. Compact neutron spectrometer 

PTB has developed the CNS based on the requirements of the JET project. It consist of a BC501A 
standard cell of small size from Bicron (0 = 2.5 cm, h = 2.5 cm) coupled to a Phillips XP2020 
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photomultiplier tube (PMT) via a light guide. The cell is assembled in a sealed aluminium housing. 
A pulsed LED is integrated into the design and used for gain correction of the PMT. A diagram of 
the detector is presented in Fig. |IJ. 
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Figure 1. A diagram of the compact neutron spectrometer, illustrating the operating principle. 

The CNS makes use of a digital acquisition system made by ENEA to satisfy JET needs Q9|, OH . 
This board allows full digitalization of the PMT signal, with a sampling rate of 200 Msamples and 
14 bit/ch resolution. The full digitalized pulses are recorded on a hard drive, which allows the 
analysis to be carried out after the experiment. The board also provides a TOF mode, using the 
PMT signal as start trigger. Neutron and y induced events can be separated using pulse-shape 
discrimination (PSD) techniques |fl 0, IH EH- 

3.1 Pulse height resolution and light output function 

The scintillator used for the CNS, BC501A (also known as NE213), is an organic compound made 
from xylene (CgZ/in) and naphthalene (C\qH%). The light output due to neutron interactions in the 
material has been studied at PTB over the past 30 years IH3T1 . The main light component comes 
from the slowing down of elastic scattered protons, but light can also be produced by the recoil of 
a -particles, electrons, or carbon nuclei. The light output is normalized to the light produced by a 1 
MeV electron (MeV ee ). 

The Monte Carlo code NRESP [11411 has been used successfully at PTB to simulate pulse height 
spectra (PHS) for a given BC501A cell geometry and neutron energy. The simulated PHS calcu- 
lated by NRESP have a shaip edge which in actual measurements is smoothed out. The simulated 
PHS therefore have to be adjusted to fit the measured PHS, and this is done using Gaussian broad- 
ening (further adjustments that are also required for a good fit are discussed in section |5T4D . The 
amount of Gaussian broadening defines the intrinsic pulse height resolution of the spectrometer. 
The pulse height function L(E) and the pulse height resolution function dL/L are characteristic of 
each individual detector. Eq.(p.ip gives a formula for dL/L in terms of fit parameters A, B and C 




(3.1) 
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The parameters can be associated with contributions that originate from different processes: the 
light transmission from the scintillator to the photocathode (A), the statistical nature of the light 
production (B), and electronic noise (C). 



4. PTB facility 

The measurements needed for the characterization of the CNS were carried out at the PTB Ion 
Accelerator Facility (PIAF) [0] (see Fig. 0). Several monoenergetic neutron beams (8, 10, 12 and 
14 MeV) and a broad neutron field (E max ~ 17 MeV, also called white field), produced using a 
cyclotron, are available. This allows a full characterization of the spectrometer in the region of 
interest (from ~ 1.5 MeV to ~ 17 MeV). 




Figure 2. The PTB experimental area, with the van de Graaff accelerator and the cyclotron 

The PTB neutron reference radiation fields used for these measurements are produced using 
D(d,n) and 9 Be(p,n) reactions. The first reaction is produced by means of a deuterium gas target. 
The beam has a small energy loss in the target (typically less than ~ 100 keV) and is fully stopped 
in the backing yielding quasi-monoenergetic neutron beams. The second reaction is produced by 
means of a thick 3 mm beryllium target, which is thick enough to fully stop the incoming protons 
yielding a neutron beam with a broad energy spectrum (white field). 

The main experimental campaign was carried out in two steps; the first involved the measure- 
ment of four monoenergetic beams (8, 10, 12 and 14 MeV), the second a measurement in the white 
field for approximately 72 hours. The selector frequency of the cyclotron was ~ 900 kHz, with a 
current of / = 900 nA on the deuterium gas target. The detector was placed at a distance of 12.29 
m from the target. The uncertainty in the distance determination is assumed to be ± 1 cm. 

In addition to the main campaign, a second campaign was carried out with neutron reference 
fields generated with the Van der Graaff accelerator. These neutron fields are produced with a 
tritiated titanium (Ti(T)) target, using the reaction 3 H(d,n) 4 He for the 14 MeV field and the reaction 
2 H(d,n) 3 He for the 2.5 MeV field. With measurements in these fields, it is possible to obtain an 
independent test of the neutron response matrix and to verify the neutron energy resolution after 
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unfolding as the neutron energies of these fields (for some particular reaction conditions) are not 
determined by the properties of the targets or the ion beams HTSIl . 

5. Method of analysis 
5.1 Time of flight method 

In the time of flight method, the kinetic energy of a particle is determined using measurements of 
the travel time over a fixed distance determined by two reference points flT2l . ITHl . 

The cyclotron provides a signal for each particle bunch that is emitted. When the acquisition 
board records a pulse from the PMT, the acquisition switches from the PMT to the signal from the 
cyclotron after a pre-defined amount of samples. A small negative pulse, caused by a difference 
in the baseline values of the two input channels of the fast filter/switch unit, indicates the switch 
between signals. The cyclotron signal is recorded and it shows up as a peak for each particle bunch 
emitted from the accelerator. Fig. |3] shows an example of this procedure. 
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Figure 3. A pulse from the PMT together with cyclotron signals recorded by the ENEA DPSD system. The 
pulse from the PMT is on the left. The TOF window starts after 120 samples (each sample is 5 ns). The 
signal polarity is inverted. 

The schematic layout of this procedure, which follows l!T7ll . is shown in Fig. |j. Note that the 
negative pulses from the PMT and reference signals are inverted by the DPSD acquisition system. 

The signals with 5 ns sampling time were interpolated to 1 ns steps and a 40% digital constant 
fraction method was used to determine the peak position for TOF evaluation, both for the pulse 
of the PMT and for the reference pulses from the cyclotron. The time difference between the 
PTM signal pulse and the first cyclotron reference pulse is the TOF value of the event. The time 
scale is based of the sampling frequency of the ADC (5 ns/sample). The time distance between 2 
cyclotron reference pulses determined in the same way is the inverse of the repetition frequency of 
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Figure 4. Principle of data acquisition for TOF determination. If an event is identified as a signal over 
threshold by the DPSD system, an output from the DPSD board is used to switch the input of a fast analog 
switch device (PTB design). This switch is delayed 500 - 600 ns to have sufficient time to sample the 
complete event (neutron, gamma, or LED pulse). After the event, the reference signal from the pickup coil 
is sampled for 2 - 3 jJ,s. Finally, the switch goes back to input 1 and the system is ready for the next event. 



the cyclotron, which is determined by comparison with a standard frequency traceable to the PTB 
atomic clock. Using this procedure, the sampling frequency of the ADC was confirmed to be 5 ns 
with a relative uncertainty of 1.5 x 10~ 4 . 

Fig. shows the time distribution of events recorded for measurements in a neutron field of 14 
MeV. For the fastest particles (i.e., photons), the TOF value is higher because the PMT pulse is used 
as the start signal. The calibration of the zero point of the TOF scale for neutrons is done using the 
TOF of photons, which travel at a fixed velocity. The time resolution achieved with this methods is 
shown by the time resolution achieved for the gamma signals (the FWHM of the gamma TOF peak 
is 2.8 ns). This proves that the sample frequency of 200 MHz that was used for the experiment 
is adequate for our application. From the uncertainty in the time determination and the assumed 
uncertainty in the distance of 1 cm, an uncertainty in the neutron energy of 200 keV at 14 MeV can 
be estimated. 

5.2 n- 7 discrimination 

The shape of the light emission from BC501A is different for signals due to neutron and gamma 
radiation, as shown on Fig. This permits the use of pulse-shape discrimination (PSD) techniques 
to distinguish between these two types of radiation. Events are discriminated using this technique 
and only neutron events are used for constructing the response matrix. 

The PSD method consists of denning short, long, and total integration gates (see Fig. §), taking 
the ratio of the short integration gate to the long integration gate, and plotting this value against 
the total integration gate, as shown in Fig. [7]. For the CNS, a value of 3/14 was selected for the 
short over long gate ratio. The PSD procedure has been implemented by the ENEA group with the 
digital board analysis software (DPSD) [fl 0, EJ. 

5.3 Satellite correction 

The reaction 9 Be(p,n) is used to obtain the white field. When a pulsed proton beam is produced, 
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Figure 5. Time-of-flight measurement of a 14 MeV monoenergetic neutron field produced at the PTB 
cyclotron (distance of 12.29 m from the target), showing (from right to left) the gamma peak, the 14 MeV 
neutron peak and the broad distribution of breakup neutron. 
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Figure 6. Pulse shapes for electrons and protons. Different integration gates are also shown. 



some of the proton bunches are not completely deflected by the internal beam pulse selector. This 
creates ghost images of the main data structure, commonly known as satellite peaks, in the pulse 
height versus TOF distribution. An illustration of this situation is given in Fig. || These events 
(in the order of 7% of the total number of events) need to be removed from the analysis because 
their TOF cannot be evaluated properly. This is achieved using a stochastic procedure in which a 
fraction of events are removed at random according to the probability that an event may belong to 
a satellite peak. 
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Figure 7. PSD identification plot. The line drawn in red separates y events (top) from neutron events 
(bottom). For this example, data from measurements made in the 12 MeV monoenergetic neutron field were 
used. 
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Figure 8. Pulse height versus TOF distribution for measurements in the white field. Only events identified 
as neutron are displayed. The red contour indicates the selection of the events with correct time evaluation. 
Events from other pulses (satellites) that overlap with this region are removed using a stochastic procedure. 

5.4 Adjustment of simulated response functions and pulse height resolution 

We use an iterative procedure to adjust the simulated response functions (i.e., the PHS that result 
from the simulation) and determine the light output function of the CNS. Each iteration consists 
of two steps. The first step is to calculate simulated response functions for neutrons with the code 
NPvESP using an estimate of the light output function (this estimate comes from known light output 
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functions for similar detectors). The second step is to adjust the simulated response functions to 
fit the measured data (using the procedure that is described in detail below), which results in a 
corrected, improved light output function that can be used for a new calculation if needed. This 
iterative procedure converges quickly (in the case of the CNS, only two iterations were needed) 
and the end result is a set of corrected simulated response functions that fit the measured data well. 

We now describe the procedure used to adjust the simulated response functions. For each neu- 
tron energy, the simulated PHS needs to be scaled in the x-axis (to account for differences in the 
gain assumed in the simulation and the gain used for the measurement) as well as in the y-axis (to 
account for differences between the neutron fluence assumed for the simulation and the neutron flu- 
ence of the measurement). It is also necessary to adjust the simulated PHS to allow for a difference 
due to a zero offset in the measurement (the zero offset is the same for all neutron energies). In 
addition to the calculated and measured neutron PHS, gamma PHS (e.g., from 207 Bi) are needed to 
determine the common PHS energy scale in equivalent electron energies. The PHS for scintillation 
detectors from photons can be calculated with the PHRESP code flU. 2()7 Bi is a gamma source 
that emits three different, well separated gamma lines. By comparison to a gamma PHS measured 
under the same experimental conditions that apply to the neutron data and a calculated gamma 
PHS, the electron equivalent energy scale and the common zero offset can be determined. This is 
achieved by including the measured and calculated gamma PHS into the analysis to determine the 
zero offset for all the PHS. A change in the neutron light output function shows up as a residual 
difference between the gamma and neutron x-axis stretching factors. The simulated PHS needs to 
be Gaussian broadened, where the Gaussian broadening is expressed by Eq. (??). 

In previous works, the fitting procedure used to adjust the simulated PHS was sequential; 
i.e., one neutron energy at the time, starting with the highest neutron energy available. In the new 
method that is applied here, all the parameters are determined simultaneously, using a maximization 
procedure based on the relative entropy calculated between measured PHS and corrected simulated 
PHS; i.e., the set of parameters that maximizes this relative entropy is the optimal set of parameters 
which gives the best fit. We used the relative entropy S in the form due Skilling [1191 EOl . 



where i m j n and i max are the minimum and maximum bins of the measured and corrected simulated 
PHS, 0J 6 and <^ s . The values of i m i n and i max are selected for each energy and correspond to the 
region of the edge on the pulse height spectrum. The maximization is carried out with respect 
to the parameters that are used to adjust the simulated spectrum. To carry out the maximization, 
we used the Broyden-Fletcher-Goldfarb-Shanno (BFGS) algorithm flTR! that is implemented in the 
GNU Scientific Library (GSL) [22D- More details of this particular application of relative entropy 
maximization can be found in [0]. 

In previous works, the adjustment of the simulated response functions was only carried out for 
those energies for which there were measurements of monoenergetic beams available. In the new 
method used here, we use the TOF information of events from the measurement in the white field 
to select only those events that correspond to a given energy (more precisely, within a small energy 
interval about the given energy). This is possible because all the required information is stored by 




(5.1) 
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the digital acquisition. A similar method has been used previously but with an analog acquisition 
system HT5I1 . 

For the characterization of the CNS, 30 neutron energies between 1.5 MeV and 17 MeV were 
selected. The energy width of each selection was chosen taking into consideration the energy 
binning that was used for the NRESP calculation and the counting statistics available. A width that 
corresponded to two of the energy bins used for the NRESP calculation was selected in order to 
have enough statistics in each of the measured PHS. The energy binning used was non- linear, it 
varied from AE 2 .5 = 23 keV at 2.5 MeV to AE14 = 66 keV at 14 MeV. 

The software written to carry out the fitting procedure provides a real-time display of the 
maximization status and, if needed, the user can easily restart the maximization in case of non- 
convergence. An example of a fit for a gamma PHS from a 207 Bi /-calibration source is shown in 
Fig. 0. 




L /MeVee 

Figure 9. The upper figure shows the measured PHS from a 201 Bi y-source (black) and the PHS simulated 
using PHRESP (red). The lower figure shows the corrected simulated PHS (red) ploted with the same 
measured PHS (black) shown in the upper figure. 

5.5 Experimentally determined response functions 

The experimental response matrix is put together from measurements made in the white field, by 
selecting data from particular neutron energy ranges. Fig. [K^ shows all neutron events which were 
analyzed (28 x 10 6 events). The energy ranges chosen correspond to the neutron energy bins used 
for the NRESP simulations flT3H . To correct for low measurement statistics, an additional step is 
carried out and smoothing is introduced. 

5.6 Unfolding of measurements of monoenergetic beams and energy resolution 

The PHS from measurements of monoenergetic neutron beams were unfolded using the code 
MAXED together with the experimental and the corrected simulated response matrices. This re- 
sulted in neutron energy spectra for the monoenergetic neutron beams. To determine the neutron 
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Figure 10. Data used for the determination of the experimental response matrix (to reduce the number of 
matrix elements plotted in the figure, a coarse binning has been used where 8 neutron energy bins and 20 
PHS bins have been combined). A total of 28 x 10 6 neutron events were evaluated. 

energy resolution, we compared the FWHM of simulations of the neutron energy beams with the 
FWHM of the unfolded spectra. We also compared to the FWHM determined using TOF. The 
results are presented in the next section. 

6. Results 

6.1 Pulse height resolution 

The procedure used to adjust the simulated spectra gives the following values for the Gaussian 
broadening parameters: A = 2.82, B = 10.0 and C = 0.6. A plot of the light output resolution (see 
eq. [Ol ) is shown in Fig. [TT]. 

The set of parameters that regulates the scaling of the simulated PHS in the x-axis for different 
neutron energies (see the discussion in subsection |5T3D is also of interest. When the light output 
function used for the NRESP simulation is adequate, this set of scaling parameters should differ 
from 1 by less than a percent. This is indeed the case after the last iteration of our procedure is 
carried out: At E n = 2.5 MeV and E n = 14 MeV, the difference from 1 is -0.90% and +0.35% 
respectively. 

6.2 Measurements in the white field 

The white neutron field provides a good test for the CNS, since it is a continuous spectrum that 
covers all neutron energies to up to ~ 17 MeV. Fig. [T2| shows unfoldings carried out using both 
experimental and corrected simulated responses, together with the result obtained using the TOF 
evaluation. 
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Figure 12. Measured spectral neutron fluence of the broad neutron field unfolded using the experimental (in 
red) and corrected simulatred (in blue) response matrices. These unfoldings are compared to the experimen- 
tal evaluation of TOF (in green). 



6.3 Mono-energetic beams 

Table [I] shows a comparison of the three different methods of estimating the mean energy and the 
FWHM of the monoenergetic neutron beams that were measured during the experimental cam- 
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Table 1. Comparison of three different methods of estimating the mean energy E n and the FWHM of the 
monoenergetic neutron beams (Gaussian fits of the neutron energy spectra with the MINUIT code were 
used to determine these parameters). The resolution (Res = FWHM /E n ) is also shown. The superscripts 
exp and sim denote the unfoldings done with the experimental and corrected simulated response matrices, 
respectively. The label TOF denotes the experimental evaluation done with time of flight. In the cases were 
TOF results were not available, the results are from the simulation code TARGET and this is indicated with 
an asterisk. 

paign. For the first set of measurements, which were carried out using the cyclotron, TOF infor- 
mation is available. The last 2 measurements were carried out with the Van de Graaff accelerator, 
and no TOF information is available. For the last two lines, we provide estimates of the mean 
energy and the FWHM derived from simulations done with the code TARGET (these values 
are indicated with an asterisk). 

In general, the unfoldings were done using a range that excluded the lower energy channels 
of the PHS but large enough to result in an unfolding that showed the complete peak. The one 
exception is the unfolding that corresponds to the values in the first line (in italics), which extends 
to much lower PHS channels so that neutrons that result from the D break-up in the d(D,n) reaction 
for neutron energies > 8 MeV also appear in the unfolded spectrum (see Figs. |T3| and 

The neutron energy spectrum of Fig. |T3| shows the main energy peak at 14 MeV and neutrons 
from the break-up process at lower energies. This illustrates clearly the good capabilities of the 
unfolding technique to reveal the whole energy distribution of the field. 

Fig shows the measured PHS together with the predicted PHS that results from folding 
the neutron energy spectrum shown in Fig. [T3] with the response functions. There is very good 
agreement between both PHS. 

Finally, a comparison was made between the neutron flux determined by the CNS and by an 
independent method (based on the PTB accelerator's reference monitor) for a monoenergetic beam 
with neutron energy E n = 14 MeV. The result derived from the PTB monitor is <P = 168600 cm" , 
while the CNS result is <P = 166200 cm~ 2 . The agreement is very good, with a difference of less 
than 1.5%. 

6.4 Comparison between neutron fields 

Fig. [T3] shows a comparison between two PHS. The first PHS, in red, is a measurement of a 
monoenergetic beam with neutron energy E n — 14 MeV. The second PHS, in black, was constructed 
from the measurement of the white field using TOF information to select only those neutrons with 



-13- 




f„ / MeV 



Figure 13. Unfolding of the 14 MeV monoenergetic field (red), showing both the main energy peak at 14 
MeV and the break-up process at lower energies. The green curve is the Gaussian fit used to extract the 
mean energy and FWHM of the peak. 




Figure 14. Measured PHS (red) compared to the predicted PHS (blue) that results from folding the neutron 
energy spectrum derived with the unfolding with the response functions. 

energies 13.985 MeV < E n < 14.015 MeV. The low statistics of the energy slice of the white field 
explains the large statistical variation observed for the PHS plotted in black. 

The excellent agreement between the two PHS provides a check of the whole procedure used 
for the analysis, in particular of the TOF evaluation that has been implemented. This comparison is 
important because the experimental response matrix is put together from PHS that are constructed 
using the energy information obtained from the TOF procedure. Also, different energy slices of the 
white field measurements were used to evaluate the pulse height resolution. 
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Figure 15. Comparison of pulse height spectra from the 14 MeV monoenergetic field (in red) and white 
field (in black). The energy selection in the white field range from 13.985 MeV to 14.015 MeV. 

7. Conclusions 

We have carried out the characterization of a CNS based on a BC501A fast-neutron detector with 
a digital acquisition system. To achieve this goal, we developed a new analysis software written in 
C++ that uses the ROOT framework. The analysis allowed us to determine two response matrices 
for the system, an experimentally determined response matrix and a response matrix based on 
simulations adjusted to fit data collected during a measurement campaign carried out at the PTB 
accelerator facility. We have evaluated both the pulse height resolution and the energy resolution 
that is achievable after unfolding the measurements with the MAXED code, and compared to the 
energy resolution determined by TOF methods. 

The results obtained for the pulse height resolution are in agreement with values that were 
obtained previously with similar systems that have larger cell sizes [0]. 

The energy resolution, as determined by the unfoldings of the measurements of monoener- 
getic beams, shows that a resolution that is comparable to the TOF resolution has been achieved. 
In particular, the resolution achieved at 2.5 and 14 MeV (see lines 2 and 6 of table [I]) satisfy the 
requirements listed in section @. These results demonstrate the good quality of the response func- 
tions and the capabilities of CNS. The energy resolution determined from the unfoldings done with 
the experimental response matrix are comparable to the energy resolution obtained with TOF, but 
not as good as that of the unfoldings done with the corrected simulated response matrix. 

The comparison between the unfoldings of the white field using the experimental and the 
corrected simulated response matrices (see Fig. [T^) show good agreement. There is structure at 
about 10 MeV in the unfolding done with the corrected simulated response matrix which is not in 
the unfolding done with the experimental response matrix or in the result derived using TOF This 
seems to be connected to the cross section libraries used for the simulation, and it is still under 
investigation at this time. 
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